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These conference proceedings cover various aspects of neutrino propagation through the 
high matter and neutrino densities near the proto-neutron star in a core-collapse supernova. 
A significant feature of this regime is the impact of neutrino-neutrino interactions. Prop- 
erties of this non-linear many-neutrino system are discussed with a particular emphasis on 
its symmetries. 

1 Introduction 

Almost all the gravitational binding energy of the pre-supernova star is converted into neu- 
trinos and antineutrinos in a core-collapse supernova, yielding a very large neutrino flux [1]. 
Consequently neutrino properties play a very important role, not only in the dynamics of a 
core-collapse supernova, but also in the r-process nucleosynthesis if supernovae are the ap- 
propriate sites [2]. Neutrinos traveling through supernovae undergo matter-enhanced neutrino 
oscillations (due to the MSW effect resulting from neutrinos interacting with the background 
electrons and positrons) much like neutrinos traveling through the Sun or the Earth [3]. How- 
ever, unlike these latter sites, in a supernova environment it is possible to have matter-enhanced 
antineutrino flavor transformations [4]. It was ascertained that, because of the large number 
of neutrinos 10^^) emitted by the proto-neutron star, background neutrinos also contribute 
to the coherent forward scattering of neutrinos in a core-collapse supernova [H [5] . Once the 
importance of the flavor-mixing non-diagonal terms coming from neutrino-neutrino interactions 
was highlighted [6], it became clear that one deals with a genuine many-body problem with 
one- and two-body interactions. 

The significance of the neutrino-neutrino interactions in core-collapse supernovae [7] and the 
possibility of the occurrence of collective effects due to those interactions were recognized early 
on [8]. In a supernova environment neutrino- neutrino interactions [6] play a crucial role both 
for neutrinos and antineutrinos [II [H El [1^ . Since such collective neutrino oscillations dominate 
the neutrino propagation much deeper than the conventional matter-induced MSW effect, it 
would impact the r-process nucleosynthesis [TTJ [T^ [T3] . There is an extensive literature on this 
subject [14] , a good starting point is several recent surveys [15j [161 HZ] • An algebraic approach 
to this problem was worked out in Ref. |18j from a many-body point of view. Such an algebraic 
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approach is helpful in exploring the hidden symmetries of the system. Hamiltonian describing 
collective neutrino oscillations possesses an SU{N) / rotation symmetry in the neutrino flavor 
space [m [191 ED] . Various collective modes, including spectral swappings or splittings arise 
from this symmetry even in the inhomogeneous or anisotropic environments |19| . One expects 
that such a complex nonlinear system may exhibit further symmetries. Indeed, several authors 
noted the presence of various conserved quantities in the collective neutrino oscillations [HJ [52] ■ 
More recently, it was shown that collective oscillations that maintain coherence can be classified 
by a number of linearly-independent functions |23| , implying that scalar products of a unique 
linear combination of the original polarization vectors are conserved. The flavor evolution of a 
dense neutrino system by taking into account both the vacuum oscillations and self interactions 
of neutrinos from a many-body perspective was considered in Ref. [24j . Using the similarity 
between the collective neutrino oscillation Hamiltonian and the BCS Hamiltonian one can show 
that, in the single angle approximation, both the full many-body picture and the commonly- 
used effective one-particle picture possess several constants of motion [53] . 

One appealing aspect of the core-collapse supernovae is that they are the laboratories where 
diverse aspects of neutrino physics come into play. Here we concentrate on the collective behav- 
ior arising from the neutrino-neutrino interactions in a supernova and omit other interesting 
topics such as the role of sterile neutrinos f26| or the neutrino magnetic moment 28; in the 
r-process nucleosynthesis, effects of the CP-violation [55], effects of turbulence and density fluc- 
tuations [30] , the role of neutrinos in shock revival |31| , and neutrino signatures of black hole 
formation [32] . 

The next chapter of this proceedings contribution describes the algebraic formulation of the 
neutrino-neutrino interactions in the many-neutrino system. Chapter 3 includes a discussion 
of the invariants of the Hamiltonian of this system. Brief concluding remarks arc included in 
Chapter 4. 

2 An Algebraic Formulation of the Neutrino-Neutrino In- 
teractions 

For simplicity, we first consider a neutrino gas with two flavors and no antineutrinos. Matter 
and flavor basis creation and annihilation operators for a neutrino with momentum p and spin 
s are related as 

ai(p, s) = cos 61 ae(p, s) — sm6 ax{p, s) (1) 
a2(p, s) = sin6' ae(p, s) + cos6 ax{p, s). (2) 
It is easy to show that the flavor isospin operators defined as 

Jp.s = o-Up, s)ax{p, s) , J",, = 4(p, s)ae(p, s) , 
Jp,s = ^ (4(P> s)ae(p, s) - 4(p, s)a^(p, s)) 

form an SU(2) algebra: 

\ ] — 2fi S P \P 1 — ±(5 

L'-'p,si •-'q.rl ^"pq"sr^p,s i L'-^p,.s' '-'q.ri ~ -^"pq"s7->^p,s ■ 

One can show that the particle mixing given in Eqs. ([!]) and ([2]) is generated by this algebra: 

U^ai{p,s)U = ae(p,s) 
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with 

U = e^p "p e' '''■^o e" ^p "p . (3) 

After subtracting a term proportional to the identity, the one-body Hamihonian including 
interactions with the electron background takes the form 

p ^ 

where one defines 

B ^ (sin26',0,-cos26'). 
The neutrino-neutrino interaction term in the Hamiltonian is 

Huv = y ^ (1 - cos dpq) Jp ■ Jq (4) 
Piq 

where ??pq is the angle between neutrino momenta p and q. In Eq. (1 — cosi?) terms follow 
from the V-A nature of the weak interactions. 
The evolution operator 

can be calculated [18] as a path integral using SU(2) coherent states: 

= exp dpz(p,t)J+(p)j|0), 10) =n4(p)|0) 

to obtain 

{z'{tf)\U\z{U)) = j V[z,z*]ci^^{iS[z,z*]) 

where 

The stationary phase approximation to the path integral 

_d _a d_\ ^, f d d d 

dt dz dz y 

yields the differential equation 



S{Z, Z)= dt + l0g(z {tf)\z{tf)). 



dtdz dzj"^^^'^*^-' UaF-a?^)^(^'^*) = « 



tzip,t) = P{p,t) - a{p,t)z{p,t) - ^*{p,t)z{p,tf (5) 

where we defined 

■l-\z{q,t)f' 



a{p, t) 



■ cos 26* + V2GFNe + V2Gf J dq{l - cosi9pJ ^i- 



and 



Z Zp J 
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Defining 



with the auxihary condition IV'eP + I'/'a-P = 1, Eq. ([5]) reduces to 



d 
'di 



A + £> - A cos 26* 
Dae + A sin 26 



-A-D 



A sin 26* 
+ A cos 20 



(6) 



where 



6m 



A = V2GFNe, 



D = V2Gf J dqil-cos^pg)[{\MQ,t)\^ ~\i'xiq,t)\'')], 

2V2Gf I dq{i - cosdp,) {Mq,tWx{q,t)) ■ 



De 

In the stationary point approximation to the full quantum mechanical problem, the test 
neutrino interacts with an "average field" representing the effect of all the other neutrinos. 
This approximation is analogous to the random phase approximation (RPA), widely used in 
many-body physics. In the RPA one can approximate product of two commuting arbitrary 
operators Oi and O2 as 

6,62 ^ Oi(cio2ie> + (eiOiie>02 - (eiOiio(^p2io, 

provided that 

{i\d,d2\0 = mi\£,){£,\d2\0. 

This approximation reduces Hui, to a one-body Hamiltonian; 
^V2Gf 



(7) 



V 



d^p d^'qR 



pq 



Mp){Mq)) + \j+{p){J-{q)) + \j-{p){J+{q)) 



The pre-exponential determinant obtained in the stationary phase approximation to the path 
integral is rather complicated and an explicit evaluation is not yet available in the literature. 
For simplicity in the discussion above we omitted antineutrinos. Antineutrinos can be included 
by introducing a second set of SU(2) algebras [18]. Similarly incorporating three flavors requires 
introduction of SU(3) algebras [Sj- Both extensions are straightforward, but tedious. 

Introducing the polarization vectors 



with the density matrix 



Px. 



P,(g)=Tr(J,((z)p) 
1 



(8) 



Pex 
Pxx 



(PO + P • fT) , 



one can show that, in RPA and including antineutrinos, the evolution equation ^ takes the 
commonly-used forms 



8 P 



and 



2p 

Sm^ 
17 



(sin 20X - cos 20z.) + \/2 Gf 



(sin 26lx - cos 29z) + V2Gf 



dq (1 - cos 'dpg) (Pq - Pq; 

dq(l-C0Sl5p,)(Pq-Pq) 



X P„ 



X P. 
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3 Invariants of the neutrino-neutrino interaction Hamil- 
tonian 



We consider the limit where the neutrino-neutrino interactions dominate and neglect the inter- 
actions with the background electrons and positrons. The Hamiltonian then becomes 



— B ■ J, + ^ (1 - COS^pq) Jp • Jq. 



2v ' V 

We further limit this discussion to the so-called single-angle approximation where all the "d are 
the Sc 
form 



the same. Defining /i = (1 — cost^) ^"^^ , t = ^t, and ujp = j^^^, the Hamiltonian takes the 



H = Y^ iOpB ■ Jp + J- J. (9) 

p 

This Hamiltonian of Eq. ^ preserves the length of each spin 

, 

as well as the total spin component in the direction of the "external magnetic field", B 

Co = S • J , \h,Co\^0 . (10) 



H, Lp 



It is possible to show that [55] the collective neutrino Hamiltonian of Eq. ^ has the following 
constants of motion: 

hp = B-Jp + 2T^^I^. (11) 

The individual neutrino spin-length discussed before in an independent invariant. However 
Co — ^p- The Hamiltonian itself is also a linear combination of these invariants. 

= ^ Wphp + y^Lp. 

p p 

The maximal value of the neutrino flavor isospin quantum number is J„,^^ ~ N/2, where N 
is total number of neutrinos. For example a state with all electron neutrinos is |i^e i/g • ■ • ) = 
I "/max >/max)/- K Is also casy to show that matter and flavor bases are connected with the 
unitary transformation of Eq. ([3]): |Jmax >/max)/ = C^^|>/max •/max)m- One has | Jmax ^max)m = 
rip^s aliP, s) |0) and | J,„^, - Jmax)™ = Yips "2(P' |0)> corresponding to the energies £'(+j„,.,^) = 
~ J2p ^^P + J^<^^ (J^--'^ + 1) and = J^p ^^p + ^^max ( ^max + l) , respectlvcly. Energy 

eigenvalues can be obtained by introducing the raising operator [25j 

Q±(^) -J— (cos^ 9J^ + sin 261 J° - sin^ OJ^ 

p ^p- i V 

Applying the operator to the state | J — J),„ yields a term proportional to Q'^{£,)\J — J)m and 
an additional term. Setting the coefficient of this latter term to zero gives the Bethe ansatz 
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equation J2p ^ -4 ^ ^1 ^^^^ needs to be satisfied if Q^{C)\J — J)m is an eigenstate. The 
most general eigenstate is 

16,6, ■•■U = Q+(ei)Q+fe)...Q+(6)k -J)m 

which corresponds to the eigenvalue 

a=l 

provided that the following Bethe ansatz equations are satisfied: 

y—i^ = -l+ y -J— . 

Using the polarization vector, Pp g = 2(Jp <;), of Eq. ([5]) one can write the Hamiltonian of 
Eq. © as 

H ^ H^^^ = ^ujpB ■ Jp + P ■ J. (12) 
p 

This Hamiltonian yields the Heisenberg equations of motion for the operators Jpi 

-^Jp = -i[Jp,i?'^^*] = (wpB + P) X Jp. (13) 
ctr 

Applying the RPA consistency conditions of Eq. (O to Eq. (|T3| one obtains the equations of 
motion in the RPA: 

— Pp = {^pB + P)x Pp. 

CLT 

Note that this is the single-angle limit of Eq. © only with neutrinos. It turns out that 
expectation value of the invariants, Eq. (|lip . of the exact many-body Hamiltonian 

Ip = 2{hp) = B-Pp+Y^ (14) 
is an invariant of the RPA Hamiltonian: 

Introduction of antineutrinos is again straightforward utilizing a second SU(2) algebra (denoted 
below with a tilde over the appropriate quantity) . The conserved quantities for the single-angle 
Hamiltonian with both neutrinos and antineutrinos 

i7 = ^ ujpB ■ Jp + ^ ujpB ■ Jp + (^J + j) • (^J + 



are 



hp^B-jp + 2y + 2y Aii^ 
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for each neutrino energy mode p and 



or each antineutrino energy mode. In the RPA these take the form 



and 



OJp — UJq U!p — U!q 



respectively. 

Recently a lot of attention was paid to the spectral splitting (or spectral swapping) phe- 
nomenon [21] [Ml [35l [Ml [32]. To explore the origin of this phenomenon we note that the 
expectation value of the invariant in Eq. (jlOp which can be written in terms of the expectation 
value of the invariants in Eq. ([T4| . 

E/p = B.P, (15) 

p 



is not conserved by the RPA Hamiltonian, Eq. ([12)) . Its conservation needs to be enforced using 



a Lagrange multiplier. Since ^ Ip is proportional to J*^, one needs to diagonalize the quantity 



p 

= E2Ap[/'tjOf7', 



p,s 

where the transforming operator is parameterized as 
with 



Zp = e^^ tan Op 



and 

cos 0. 



This operator transforms matter-basis creation and annihilation operators into quasi-particle 
creation and annihilation operators: 

ai(p, s) = C/'^ai(p, s) ?7' = cos 0p fli (p, s) — e'* sin 6'p 02 (p, s) 
a2(p, s) = t/'^a2(p, s)f/' = e~"' sin 6'p ai(p, s) -I- cos 6'p 02 (p, s) 
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so that we obtain a diagonal Hamiltonian: 

7?"^* + c^, J" - ^ Ap (al(p,s)ai(p,s)-4(p,s)a2(p,s) 



Let us assume that initially (lim /i — > oo) there are more Ve^s and all neutrinos are in flavor 
eigenstates. We then have 



lim cos 9p = lim 



po ^ 
— cos 29 
\P\ 



cos 6. 



i.e., the diagonalizing transformation of Eq. (jl6p reduces into the neutrino mixing 
tion of Eq. (|3]) and the total Hamiltonian of Eq. (|T6)) is diagonalized by the flavor 

ai(p,.s) U^ai{p,s)U = ae(p,s). 

After neutrinos propagate to a region with very low neutrino density (/z — > 0) one 



transforma- 
eigenstates: 



gets 



cos 9^ 



\uJc - UJp\ 



1 Up < LOf. 
> CJc 



yielding 



ai(p,s) = C/Ta^(-p^s)[/ ^ J ^^i^H, ; p 



-a2(p,s) Wp > Wc ' 

i.e. neutrinos with lOp < Uc and ujp > uic evolve into different mass eigenstates. In 
was shown that such an evolution leads to spectral splits. 



Ref. [2T] it 



4 Concluding Remarks 

Neutrino propagation through the dense media in the core-collapse supernovae probes many 
interesting collective effects. Because of the neutrino-neutrino interactions, this many-body sys- 
tem is intrinsically non- linear, it can be linearized only in certain cases [35] ■ We examined this 
many-neutrino system both from the exact many-body perspective and from the point of view 
of an effective one-body description formulated with the application of the RPA method. To 
achieve this goal we exploited mathematical similarities between the neutrino-neutrino interac- 
tion Hamiltonian and the BCS pairing Hamiltonian. (Indeed, the N-mode collective oscillations 
of the neutrinos are related to the m-spin solutions of the BCS model [39]). In the limit of 
the single angle approximation, both the many-body and the RPA pictures possess many con- 
stants of motion manifesting the existence of associated dynamical symmetries in the system. 
Judicious use of these invariants could certainly help numerical calculations |40| . 

We treated the two-flavor neutrino-neutrino interaction in the single-angle approximation, 
and provided an interpretation of the critical energy in the spectral swap/split phenomenon as 
the Lagrange multiplier of the number conservation constraint. Recent numerical work with 
three flavors in the multi-angle approximation uncovers significant differences between single- 
and multi-angle formulations j41j . In particular, multi-angle formulation is found to reduce the 
adiabaticity of flavor evolution in the normal neutrino mass hierarchy, resulting in lower swap 
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energies. Thus it seems that single-angle approximation seems to be sufficient in some cases, 
but inadequate in other situations. 

Other questions remain regarding the many-body behavior of the neutrino system. For 
example, in the calculations so far neutrinos are assumed to be emitted half-isotropically (only 
outward- moving modes are occupied with backward- moving modes being empty). However, 
recent realistic supernova simulations suggest that neutrino angular distributions arc not half- 
isotropic [42] . Flavor-dependent angular distributions may lead to multi-angle instabilities 
j43U44| . Future work could uncover even more interesting features of this many- neutrino system. 
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